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ABSTRACT 

A  dynamic  simulation  model  which  has  been  specifically 
designed  to  determine  the  hydrologic  changes  resulting  from 
timber  harvesting  and  correlary  models  which  simulate  timber 
yields  are  described.  Emphasis  is  placed  on  the  "planning  unit" 
which  is  defined  by  environmental  characteristics,  including 
combinations  of  slope,  aspect,  elevation,  and  forest  cover.  The 
models  are  intended  for  use  on  subalpine  watersheds  where  the 
primary  source  of  streamflow  is  from  melting  snow.  The 
hydrologic  model  simulates  winter  snow  accumulation,  the  short- 
and  longwave  radiation  balance,  snowpack  condition,  snowmelt, 
and  subsequent  runoff  in  time  and  space.  The  timber  models 
simulate  projected  timber  yields  in  response  to  changes  in  cultural 
treatments  and/or  variations  in  original  stand  and  site  con- 
ditions. 

The  models  are  capable  of  simulating  a  broad  array  of  timber 
harvesting  alternatives.  Hydrologic  changes  and  timber  yields 
can  be  determined  for  intervals  of  time  which  can  vary  from  a  few 
years  to  the  rotation  age  of  subalpine  forests  (120  years  and 
longer).  In  the  hydrologic  model,  this  is  accomplished  by  means  of 
time  trend  functions  which  compute  changes  in  evapotranspira- 
tion,  soil  water,  forest  cover  density,  reflectivity,  interception,  and 
snow  redistribution  as  the  forest  stands  respond  to  management. 

The  models  have  been  used  to  simulate  the  effects  of  forest  and 
watershed  management  on  several  representative  drainage 
basins  in  the  Rocky  Mountain  Region  of  the  United  States. 
Projected  hydrologic  changes  and  growth  and  yield  subsequent  to 
timber  harvesting  in  lodgepole  pine  and  spruce-fir  are  described  in 
this  report. 
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Introduction 

Watershed  management  research  during 
the  past  50  years  has  shown  that  subalpine 
forests  exert  a  significant  effect  on  water 
yields.  Hence,  man-caused  changes  in  the 
forest  environment  can  be  expected  to  affect 
the  water  resource.  When  timber  is  harvested, 
the  magnitude  of  the  resulting  hydrologic 
change  is  highly  sensitive  to  the  pattern  in 
which  a  given  volume  of  wood  is  removed. 

Because  water  from  snowmelt  is  a 
primary  resource  in  the  Rocky  Mountain 
Region  of  the  United  States,  the  need  for  a 
planning  tool  to  evaluate  the  potential 
hydrologic  effects  of  various  timber 
harvesting  strategies  is  obvious  in  view  of  the 
U.S.  Forest  Service  goal  of  sound  multi- 
resource  management. 

Some  progress  has  been  made  in  the 
development  of  dynamic  simulation  models 
which  predict  the  short-term  effects  of  timber 
harvesting  on  snowmelt  and  water  jdeld 
(Leaf  and  Brink  1972,  1973a,  1973b).  This 
work  has  recently  been  expanded  to  deter- 
mine the  long-term  interactions  between  the 
water  and  timber  resources  with  regard  to 
initial  partial  cutting  in  old-growth  subalpine 
forests,  and  subsequent  management  of  these 
stands. 

The  system  described  in  this  report 
utilizes  output  from  a  water  balance  model 
(Leaf  and  Brink  1973b,  1975)  to  simulate  both 
the  immediate  and  long-term  effects  of  forest 
and  watershed  management,  in  areas  where 
runoff  is  derived  primarily  from  melting 
snow.  Considerable  flexibility  is  provided  for 
simulating  alternative  silvicultural  systems. 
Moreover,  yield  tables  can  be  produced  which 
show  how  projected  timber  volumes  will  vary 
in  response  to  various  timber  management 
alternatives  (Myers  1971). 

The  hydrologic  model  is  described  first, 
followed  by  examples  of  projected  hydrologic 


changes  from  timber  harvesting  on  two 
watersheds  in  Colorado  and  Wyoming.  A 
discussion  of  the  timber  models,  which 
simulate  growth  and  yield  once  old-growth 
subalpine  forests  are  converted  to  managed 
stands  is  also  presented.  The  models  simulate 
stand  growth  and  response  to  intermediate 
cuttings  from  the  regeneration  period  to  final 
harvest  for  timber  production.  Finally,  ex- 
amples are  given  which  illustrate  how  the 
models  can  be  used  simultaneously  to  provide 
the  manager  with  multi-resource  response 
data  for  timber  and  water. 


Hydrologic  Model 
Model  Configuration 

The  analytical  framework  of  the  system  is 
based  on  a  "planning  unit"  which  is  defined 
by  environmental  characteristics  including 
combinations  of  slope,  aspect,  elevation,  and 
the  species,  form,  and  structure  of  the  forest 
cover.  The  model  is  designed  to  simulate  the 
hydrologic  effects  of  timber  harvesting  in 
order  to  develop  management  strategies  for 
planning  intervals  which  can  vary  from  a 
few  years  to  the  rotation  age  of  subalpine 
forests  (120  years  and  longer). 

Management  strategies  may  subdivide  a 
given  planning  unit  into  as  many  as  eight 
distinct  areas  or  "response  units,"  which 
may  be  managed  independently  at  varying 
points  in  time  during  the  planning  interval. 
Provision  is  also  made  so  that  different 
cutting  practices  may  be  imposed  on  the 
response  units,  and  finally,  any  number  of 
cuttings  may  be  made  on  a  given  response 
unit  at  specified  years  during  the  planning 
interval. 


Hydrologic  integrity  is  maintained  as 
management  strategies  are  formulated,  since 
all  interactions  between  the  various  response 
units  are  accounted  for  in  time  and  space. 
The  interactive  effects  of  a  new  decision  on 
ones  previously  implemented  are  simulated, 
as  are  the  effects  of  time,  as  demonstrated 
through  reforestation.  Moreover  the  overall 
hydrologic  effects  resulting  from  each 
management  decision  on  the  planning  unit 
are  projected  to  the  end  of  the  planning 
interval  as  though  it  were  the  final  decision 
in  the  strategy.  Thus,  the  singular  effects  of 
each  decision  can  be  evaluated. 

A  core  model  simulates  the  water 
balance  on  a  daily  basis.  On  those  areas 
where  forest  cover  has  been  removed,  the 
parameters  which  define  soil  water 
availability,  forest  cover  density,  reflectivity, 
interception,  and  snow  distribution  are  ad- 
justed on  an  annual  basis  by  means  of  time- 
trend  relationships  described  later  in  this 
report. 

Because  climatological  observations  are 
rarely  available  for  the  long  periods  of  time 
simulated,  the  system  has  the  capability  to 
extend  a  sample  data  base  by  a  randomized 
selection  of  water  years  until  the  desired 
planning  interval  is  completed. 


Theory 

The  discussion  which  follows  is 
necessarily  brief  and  intended  to  give  the 
reader  a  general  idea  of  the  scope  of  the 
model.  Detailed  flow  chart  descriptions  and 
more  complete  hydrologic  theory  have  been 
pubHshed  by  Leaf  and  Brink  (1973a,  1973b, 
1975).  The  reader  is  referred  to  these 
publications  for  a  more  comprehensive  treat- 
ment of  the  concepts  which  follow.  At  the 
present  time,  the  model  is  being  used  in 
selected  areas  throughout  the  Rocky  Moun- 
tain Region.  A  flow  chart  of  the  core  model  in 
the  system  is  shown  in  figure  1.  Figure  2 
shows  how  the  core  model  is  used  in  ex- 
ecuting alternative  management  strategies. 

Snowmelt.  —  Previous  work  in  high- 
elevation  coniferous  forests  has  shown  that 
radiation  is  the  major  source  of  energy  for 
snowmelt.  Accordingly,  short-  and  longwave 
radiation  represent  the  energy  components 
available  for  snowmelt.  Shortwave  radiation 
to  the  snow  or  ground  surface  beneath  the 
forest  canopy  is  controlled  by  a  transmissivi- 
ty  coefficient  to  be  discussed  later,  and  which 


varies  according  to  forest  cover 
characteristics.  The  incident  shortwave 
radiation  as  measured  on  a  horizontal  sur- 
face is  adjusted  according  to  the  slope  and 
aspect  of  each  hydrologic  response  unit. 
Longwave  radiation  is  computed  by  the 
Stefan-Boltzmann  equation. 
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Figure  1 
model. 


—  Flow  chart  showing  system  core 


Snowpack  reflectivity  is  varied  accord- 
ing to  precipitation  form,  air  temperature, 
and  the  energy  balance.  During  the  winter 
months,  temperatures  within  the  snow  cover 
are  simulated  using  unsteady  heat  flow 
theory.  The  snowpack  will  yield  melt  water 
only  when  it  has  become  isothermal  at  0°C, 
and  its  free  water  holding  capacity  is  reach- 
ed. 

The  generalized  equation  for  snowmelt 
can  be  written  as: 

M  -  S^d.O  -  R)T  +  C^(aT^4  _  „Tg^) 

+  (1.0  -C^XaaT^^  -aTg^)  [i] 
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where 

M  =  Daily  snowmelt  in  calories  per 
cm^, 

S\Y  =  Incoming  solar  radiation  in  ly  per 
day; 

R  =  Reflectivity  of  snowpack,  express- 
ed as  a  decimal; 

T  =  Shortwave  radiation  transmissiv- 
ity  coefficient,  expressed  as  a  deci- 
mal; 

C(j  =  Forest  cover  density,  expressed  as 
a  decimal; 

a      =  Stefan-Boltzmann  Constant; 

Ta    =  Air  temperature  in  °K; 

Tg    =  Snow  surface  temperature  in  °K; 

a  =  Coefficient  for  computing  sky  radi- 
ation, expressed  as  a  decimal.  On 
clear  days,  a  =  1.00. 


Evapotranspiration.  —  A  "potential" 
evapotranspiration  function  was  developed 
for  the  model  based  on  the  empirical  Hamon 
equation  (Hamon  1961),  which  requires 
latitude,  converted  to  day  length;  and  mean 
monthly  temperature,  converted  to  satura- 
tion vapor  density.  The  coefficient,  C,  in 
Hamon's  equation  was  modified  in  order  to 
obtain  an  expression  for  potential 
evapotranspiration,  Eg  under  "unlimited" 
solar  input,  assumed  herein  as  potential 
radiation.  The  evapotranspiration  computed 
by  this  expression  is  reduced  in  proportion  to 
the  radiation  actually  received  each  day 
according  to  the  expression: 


E, 


(C'D2p  )(?^) 
t  p 


[2] 
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Figure  2.  —  Flow  chart  showing  how  core 
model  is  used  to  execute  alternative 
management  strategies 


where 

C  =  the  modified  coefficient  defined 
above, 

D  =  possible  sunshine  in  units  of  12 
hours, 

p^  =  saturated  water  vapor  density  (ab- 
solute humidity)  at  the  daily  mean 
temperature  in  grams  per  cubic 
meter, 

SW  =  daily  shortwave  radiation  in  lang- 
leys,  and 

P  =  potential  shortwave  radiation  for 
the  day  as  computed  by  Frank  and 
Lee  (1966). 


The  adjusted  evapotranspiration  is  then 
redefined  according  to  its  source,  which  can 
include:  (1)  evaporation  from  snow  in- 
tercepted by  the  forest  canopy;  (2)  evapora- 
tion from  the  snowpack  surface;  and  (3) 
evapotranspiration  during  the  growing 
season. 

Input  to  the  watershed  system  is  derived 
from  snowmelt  and  rainfall.  Once 
evapotranspiration  requirements  have  been 
satisfied,  any  remaining  input  is  used  to 
satisfy  soil  water  recharge  requirements. 
When  field  capacity  is  reached,  the  residual 
input  becomes  water  available  for  streamflow 
(generated  runoff). 

With  regard  to  evapotranspiration,  it  was 
assumed  that  water  use  by  old-growth  forest 
during  the  growing  season  proceeds  at  rates 
limited  only  by  available  energy  until  the  soil 
water  is  depleted  to  50  percent  of  the  max- 
imum "available"  for  transpiration  (field 
capacity  index).  Thereafter,  transpiration  is 
decreased  in  proportion  to  the  amount  of  soil 
water  below  one-half  of  the  field  capacity 
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index.  In  open  or  cutover  areas,  it  was  reason- 
ed that  the  absence  of  dense  vegetation  and  a 
shallow  rooting  depth  enables  evapotrans- 
piration  to  proceed  at  maximum  rates  only 
when  the  soil  mantle  is  completely  recharged. 
Thereafter,  evapotranspiration  is  linearly 
decreased  to  zero  at  an  assumed  three-fourths 
of  the  field  capacity  index.  These 
relationships  are  shown  graphically  in  figure 
3a.  As  forest  vegetation  reoccupies  cutover 
areas,  and  consumptive  use  is  increased,  the 
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Figure  3.  —  Evapotranspiration  as  a  function 
of  available  soil  water  for:  (a)  old-growth 
forest  and  open  conditions,  and  (b)  old- 
growth  forest  and  some  intermediate  forest 
cover  condition  several  years  after  timber 
harvesting. 


relationship  in  figure  3a  changes  until  ul- 
timately, as  the  forest  cover  is  reestablished, 
it  approaches  that  of  the  old-growth  forest 
curve.  It  is  this  phenomenon  which  is 
primarily  responsible  for  diminishing  water 
yield  increases  over  time  following  timber 
harvest.  The  rate  at  which  this  transition 
takes  place  depends  upon  forest  species, 
climate,  stand  conditions,  and  the  objectives 
of  management. 

A  general  expression  for  the  relations 
shown  in  figure  3  can  be  written  as  follows: 


6  =  I 

P>  r, 
_  1 


A\P-T  +1 


<  P  < 


P  <   T  - 


e  =  0 


[3] 


where 


IS 


the 


A 


E 

the    ratio,        a  . 

evapotranspiration  rate  adjusted 
for  available  soil  water,  and  Eg  is 
computed  in  this  model  by  a  modi- 
fied version  of  the  Hamon  equation, 
the  available  soil  water  at  any  time 
during  a  given  water  year.  0  <  /? 
<  M  where  M  is  the  "field  capac- 
ity index," 

the  critical  point  at  which  avail- 
able soil  water  begins  to  limit  evap- 
otranspiration. M/2  <  r  <  M, 
and 

the  slope  of  the  relationship  be- 
tween E^/Eg  =  0  and  1. 


It  appears  from  the  Fool  Creek  watershed 
study  in  central  Colorado  that  hydrologic 
changes  resulting  from  timber  harvest  in  the 
subalpine  zone  persist  for  many  years.  The 
Fool  Creek  study  showed  that  water  yield 
increases  did  not  decrease  significantly  more 
than  16  years  after  treatment  (Hoover  and 
Leaf  1967).  These  results  and  results  from 
timber  management  research  were  used  to 
develop  the  time-trend  relationships  discuss- 
ed below.  As  seen  below,  the  procedure  used  in 
formulating  each  time-trend  relationship 
was  to:  (a)  establish  plateaus,  and  maximum 
and  minimum  values  for  each  hydrologic 
variable;  (b)  establish  critical  times  at  which 
a  transition  begins  to  occur;  and  (c)  assume  a 
functional  relationship  which  determines  all 
intermediate  values. 
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It  should  be  emphasized  that,  due  to  a 
lack  of  understanding  of  long-term 
hydrologic  phenomena,  the  time-trend 
equations  are  not  exact  in  any  intrinsic  or 
mathematical  way  (Forrester  1961).  They 
should  be  considered  only  as  relationships 
which  represent  rational  estimates  of  how  the 
most  significant  processes  vary  over  a  long 
period  of  time.  These  time-trend  relationships 
are  logically  plausible,  but  additional 
research  is  needed  before  more  precice 
equations  can  be  developed. 

Soil  Water.  —  The  critical  point  at 
which  available  soil  water  begins  to  limit 
evapotranspiration  (r),  was  assumed  to  vary 
with  time  and  forest  tree  species.  These 
relationships  were  expressed  as 


r  =  Me-k(t  -  tc,) 


r  =  M,  t  <  t„ 

T  =  M/2,  t  >  t^ 


[4] 
where 

k  =  an  index  of  the  rate  of  decline  of  r 
t^j  =  the  time  at  which  available  soil 
water  begins  to  limit  evapotrans- 
piration in  years,  and 
tj.  =  the  time  at  which  the  hydrologic 
effect  of  timber  harvesting  be- 
comes insignificant. 

The  parameters,  k  and  tcj,  vary  according  to 
tree  species.  The  assumed  relationship 
between  A  and  r  is  given  by 


M2 


[5] 


Substituting  equations  [4]  and  [5]  into  equa- 
tion [3]  yields 


=  4e 


■k(t  -  t^  ) 


Ci '  IS/M  -  e 


-k(t  -  t 


+  1 


defined  as  "canopy"  or  "crown"  closure,  but 
rather  as  a  tree  parameter  which  integrates 
the  net  effects  of  the  overstory  on  the 
transmission  of  solar  radiation  to  the  forest 
floor.  Forest  cover  density  varies  according  to 
crown  closure,  the  vertical  foliage  distribu- 
tion, species,  season,  and  stocking 
(Reifsnyder  and  Lull  1965).  Empirical 
relationships  between  various  timber  stand 
variables  and  percent  radiation  beneath  the 
forest  canopy  (transmissivity  coefficient) 
have  been  derived  for  the  three  major  sub- 
alpine  tree  species  in  the  process  of  model 
calibration  and  from  solar  radiation 
transmission  studies  in  central  Colorado.  The 
resulting  equation  from  this  work  is  given  by 


T  =  0.19C 


dmx 


0.6 


[7] 


[6] 


where 

T  =  the  transmissivity  of  the  forest  can- 
opy expressed  as  a  decimal  fraction 
of  the  amount  of  solar  radiation  avail- 
able above  the  forest  canopy,  and, 

C,  =  the  natural  old-growth  forest  cover 
dmx  density  expressed  as  a  decimal. 


Combinations  of  Cdmx  and  T  for  the  three 
major  subalpine  forest  species  are  given  in 
Leaf  and  Brink  (1975). 

As  trees  reoccupy  cutover  areas,  forest 
cover  density  (Cd)  increases  with  time  until  it 
reaches  a  maximum  value.  Research  has 
shown  that  the  rate  at  which  forest  cover 
density  reaches  this  plateau  depends  on 
environmental  conditions,  stocking  levels, 
and  species.  In  subalpine  coniferous  forests 
in  the  Rocky  Mountains,  it  can  vary  from  30 
to  more  than  80  years  (Alexander  1974). 
Accordingly,  C^  was  assumed  to  vary  as  a 
function  of  time  according  to  the  following 
equation: 


which  is  a  general  equation  for  0  as  a  function 
of  forest  cover  type,  field  capacity  index,  and 
time. 

Forest  Cover  Density.  —  Forest  cover 
density  plays  an  important  role  in  the  simula- 
tion model.  It  is  the  major  descriptive 
parameter  of  the  form,  structure,  and 
arrangement  of  forest  stands,  and  therefore 
controls  the  energy  balance,  interception, 
and  evapotranspiration.  This  parameter  is 
also  related  to  basal  area  which  is  pertinent  to 
timber  production,  as  discussed  later.  Forest 
cover  density  as  used  in  the  model  is  not 


'dmx 

<j)2 


t  -  t. 


t„  <  t  <  $ 


[8] 


where 

Cd   =  intermediate  forest  cover  density 

expressed  as  a  decimal 
^     =  the  time  in  years  from  tc^  at  which 

maximum    forest    cover  density 

(Cdmx)  is  reached,  and 


t, 


C2 


=  critical  time  at  which  regeneration 
is  sufficient  to  reestablish  the 
stand.  When  t  <      ,      =  0. 


5 


Reflectivity.  —  Studies  of  the  energy 
balance  and  associated  vapor  loss  indicate 
that  the  major  variations  with  regard  to 
latent  heat  flux  are  associated  with  reflectivi- 
ty (Baumgartner  1967).  Accordingly,  a 
relationship  between  reflectivity  and  forest 
cover  density  was  derived  as  follows: 


Rf  =  R^-^exp 


'dmx 


(t  -  t^ 


$2 


[9] 


where 


R. 


=  the  reflectivity  of  the  forest  stand, 
"^fo  ~        reflectivity  of  a  forest  opening 
(assumed    herein    as    0.5).  When 
t  <  tc2,  Rf  =  0.5,  and 
-    =  1.609  Cd^ix^. 

Equation  [9]  is  used  to  adjust  equa- 
tion [6]  for  net  available  energy. 
Thus: 


ea  -  Rf) 


[10] 


where  Ea'  is  the  actual  evapotranspiration 
adjusted  for  both  available  soil-water  and 
energy. 

Interception.  —  In  the  interception 
portion  of  the  model,  it  is  assumed  that: 

1.  The  amount  of  snow  intercepted 
varies  according  to  forest  cover  densi- 
ty, Cd; 

2.  The  intercepted  snow  rests  on  the 
canopy  for  only  1  day  following  the 
day  of  the  snow  event  due  to  turbulent 
winds  which  remove  snow  from  the 
crowns;  and 

3.  The  residual  intercepted  snow  which 
is  not  vaporized  after  that  period  of 
time  is  added  to  the  snowpack. 

The  second  and  third  assumptions  are 
based  on  field  studies  which  indicate  that 
snowfall  is  strongly  influenced  by  wind  in- 
teracting with  the  forest  and  local 
topography  (Hoover  and  Leaf  1967,  Hoover 
19;69). 

Evaporation  from  the  snow  surface  and 
from  snow  intercepted  by  the  forest  canopy  is 
computed  by  the  following  rational 
relationships  (Leaf  and  Brink  1975): 


V   =  (1  -  C  .)E^ 
s     ^        d'  s 


[11] 


V   =  -!-E 
c      Cd  s 


C, 


^dmx 


[12] 


where 


E. 


=  evaporation  from  the  snow  surface 
=  intercepted  snow  evaporation, 
=  intermediate  forest  cover  density  as 
=  defined  in  equation  [8];  and 
=  potential     evapotranspiration  as 
defined  in  equation  [2]. 


When  C  ,  > 
d 


C 


dmx 


and  snow  rests  on 


the  canopy,  evaporation  is  computed  by 
equation  [12],  whereas  during  conditions 
when  the  canopy  is  free  of  snow,  evapora- 
tion takes  place  according  to  equation  [11] 

C, 


However,  when  0  <  C 


Mmx 


and  snow 


d~  2 

rests  on  the  canopy,  both  equations  [11] 
and  [12]  are  used  as  follows: 


E, 


1  - 


2C, 


^dmx 


C 


dmx 


1  -  C, 


[13] 


where 
V,  = 


combined  evaporation  from  snow 
surface  and  intercepted  snow  in 
cutover  areas. 


Equation  [13]  more  realistically  re- 
presents the  evaporation  from  cutover  areas 
which  are  not  completely  occupied  by  trees. 
Equation  [13]  applies  only  when  Cd>0 
When  Cd  =  0,  Vt  =  V_ 
equation  [8]  into  equation  [13] 
following  relationship  is  obtained 


g.  By  substituting 


the 


Vt  =  Es 


+  1 


^dmx 


1- 


c 


dmx 

<J)2 


[14] 


which  expresses  as  a  function  of  Cdmx 
and  time. 

Snow  Redistribution.  —  Redistribu- 
tion of  snow  as  a  result  of  patch-cutting  is  a 
significant  factor  influencing  runoff. 
Moreover,  in  the  lodgepole  pine  type  in 
Colorado,  this  phenomenon  is  not  greatly 
diminished  more  than  30  years  after  timber 
harvest  in  spite  of  regrowth  of  trees  and 
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associated  increase  in  forest  cover  density 
(Hoover  and  Leaf  1967,  Hoover  1969).  It  is 
believed  that  changes  in  natural  snow  ac- 
cumulation patterns  produced  by  timber 
harvest  will  persist  until  the  new  crop  of  trees 
approaches  the  height  of  the  remaining 
virgin  forest.  Moreover,  optimum  redistribu- 
tion of  snow  results  when  old-growth  sub- 
alpine  forests  are  (a)  harvested  in  small 
patches  less  than  8  tree-heights  in  diameter; 
(b)  protected  from  wind;  and  (c)  interspersed 
so  that  they  are  5  to  8  tree-heights  apart.  More 
snow  is  deposited  in  the  openings,  and  less 
snow  accumulates  in  the  uncut  forest  so  that 
total  snow  on  headwater  basins  is  not 
significantly  increased.  Accordingly,  the 
following  relationships  were  developed  for 
simulating  snow  redistribution  effects  with 
time  and  the  three  primary  tree  species: 


where 


P  =  1,  t  >  t. 


[15] 


p  =  snow  redistribution  factor  in  the 
cutover  area  which  varies  accord- 
ing to  the  silvicultural  system  used. 

Pfnx  ~  redistribution  factor  immedi- 
ately after  timber  harvesting. 

kj     =  an  index  of  the  rate  of  decline  of  p, 

tcg  =  the  time  at  which  forest  regrowth 
begins  to  reduce  snow  redistribu- 
tion in  years,  and 

tj.^  =  the  time  at  which  forest  regrowth 
causes  snow  redistribution  to  be- 
come insignificant. 

The  parameters,  k,  tr^  and  tc^  vary 
according  to  tree  type.  When  t  <  tc;,,  no 
adjustments  are  made  in  the  redistribution, 
since  field  studies  in  the  Rocky  Mountain 
Region  indicate  that  a  correction  is  not 
warranted  for  several  years  after  harvest 
cutting. 

It  should  be  emphasized  that  redistribu- 
tion theory  is  valid  only  when  timber  is 
harvested  in  small  patches  (5  to  8  tree- 
heights  in  diameter)  which  occupy  less  than 
50  percent  of  a  given  planning  unit.  An 
optimum  redistribution  factor  is  ap- 
proximately 1.30,  which  corresponds  to  5  to  8 
H  patches  which  occupy  40  percent  of  the 
planning  unit.  In  this  situation,  the  snow- 
pack  is  increased  30  percent  in  the  openings 
and  decreased  20  percent  in  the  uncut  forest. 


When  openings  are  larger,  snow  is  scoured 
from  the  center,  whereas  smaller  openings 
also  do  not  trap  snow  efficiently. 

Individual-tree  Selection  Cutting. — 
Selection  cutting  in  the  model  corresponds  to 
a  reduction  of  the  forest  cover  density  (C(j). 
The  degree  that  Cd  is  reduced  depends  on 
characteristics  of  the  stand  and  the  volume  of 
timber  removed.  In  old-growth  stands,  if  Cd 
is  reduced  by  50  percent  or  less  from  Cdmx>  it 
is  assumed  that  forest  canopy  density  does 
not  increase  subsequent  to  harvest  cutting. 
However,  if  Cd  is  reduced  more  than  50 
percent  from  Cdmx.  equation  [8]  is  used  to 
simulate  redevelopment  of  the  canopy  with 
time.  Solving  equation  [8]  for  time  yields: 


*2C 


C 


dmx 


1/2 


+  t 


C2 


[16] 


If  the  degree  to  which  thinning  reduces  Cdmx 
is  given  by  rj,  then  Cd  is  given  by 

Cd  =  Cdmxd  -  ") 

Hence,  equation  [16]  can  be  written  as: 


[d  -  v)] 


1/2 


+  t. 


C2 


[17] 


where 

t^  =  the  time  required  to  reach  the 
reduced  forest  cover  density  as  if 
the  stand  were  initially  patch-cut, 
and 

rj  =  the  degree  that  Cd  is  reduced  from 
Cdmx  expressed  as  a  decimal. 

All  of  the  time  trend  relationships  are  then 
initialized  at  t^  in  order  to  simulate  the 
hydrologic  effects  of  selection  cutting. 


Applications 

Field  Studies.  —  Watershed  studies  in 
the  Central  Rocky  Mountains  show  that 
timber  harvesting  significantly  affects  the 
hydrologic  system.  For  example,  on  the  714- 
acre  Fool  Creek  watershed  where  39  percent 
of  the  area  was  clearcut  in  strips  1  to  6  tree- 
heights  wide  (fig.  4),  snow  accumulation, 
melt,  and  subsequent  water  yield  were  all 
affected.  Hoover  and  Leaf  (1967)  report  that 
total  snow  storage  on  Fool  Creek  did  not 
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increase  after  harvest  cutting.  Strip  cutting 
caused  more  snow  to  accumulate  in  the 
openings,  however,  and  less  in  the  uncut 
forest.  When  regressed  against  a  1,984-acre 
control  watershed,  it  was  determined  that  the 
average  annual  runoff  increased  more  than  3 
inches  after  treatment  (fig.  5).  Seasonal  peak 
flows  were  not  increased,  nor  were  summer 
recession  flows  diminished  (Leaf  and  Brink 
1972).  Timber  harvesting  caused  higher 
snowmelt  rates  in  early  spring  and  more 
efficient  water  yield. 


Figure  4.  —  Fool  Creek  experimental  water- 
shed, Fraser  Experimental  Forest.  Control 
watershed  is  to  the  right  of  Fool  Creek. 


Model    Studies    of   Snowmelt.  — 

Dynamic  hydrologic  models  are  useful  tools 
for  quantifying  the  effects  of  watershed 
changes  on  runoff.  We  have  used  this 
procedure  to  study  the  effects  of  hypothetical 
watershed  management  practices  on  un- 
disturbed watersheds  in  the  Rocky  Mountain 
Region  using  our  best  information  from  field 
studies  and  the  model  described  above. 


Short-Term  Hydrologic  Impacts  of 
Timber  Harvesting 

In  simulating  a  hypothetical  timber 
harvest  on  the  667-acre  Deadhorse  Creek 
watershed  in  central  Colorado,  the  snowmelt 
portion  of  the  model  has  produced  results 
similar  to  those  observed  from  the  Fool  Creek 
experiment.  Elevations  on  Deadhorse  Creek 
vary  from  9,450  feet  msl  to  11,600  feet  msl. 
Soils  are  derived  from  gneiss  and  schist 
rocks;  the  forest  cover  is  old-growth  lodgepole 
pine  and  spruce-fir. 

Leaf  and  Brink  (1972)  assumed  that  40 
percent  of  the  watershed  area  was  uniformly 
patch-cut  in  openings  5  to  8  tree-heights  in 
diameter.  Because  field  studies  have  shown 
that  total  snow  storage  is  not  changed  follow- 
ing harvesting,  the  snowpack  was  increased 
30  percent  in  the  openings  and  decreased  20 
percent  in  the  uncut  forest.  In  addition  to 
redistributing  the  snowpack  to  represent  the 
harvesting  system,  the  forest  canopy  density 
parameter  (C(Jmx)  reduced  to  zero  on  40 
percent  of  the  area  in  each  of  10  hydrologic 
subunits. 

Results  obtained  through  manipulating 
the  input  and  forest  cover  parameters  in  the 
calibrated  model  indicated  that  patch-cutting 
small  openings  in  mature  lodgepole  pine  and 
spruce-fir  forest  results  in  increased 
snowmelt  early  in  the  melt  season  with 
diminished  snowmelt  later.  Although  timber 
cutting  affected  the  timing  of  snowmelt,  it 
apparently  did  not  significantly  change  the 
duration  of  the  snowmelt  season.  Under  com- 
parable conditions,  snowmelt  began  a  few 
days  earlier  in  small  openings,  but  in  both  the 
natural  forest  and  cutover  areas,  the  last 
snow  melted  out  at  about  the  same  time. 
Because  melt  rates  in  openings  were  higher 
early  in  the  snowmelt  season,  peak 
streamflow  would  not  increase  appreciably,  if 
at  all,  under  the  assumed  timber  harvesting 
alternative.  Figure  6a  summarizes  the 
predicted  change  in  snowmelt  input  resulting 
from  this  practice  for  the  1947-71  period  of 
record. 

In  addition  to  redistributing  the  snow- 
pack and  accelerating  snowmelt  runoff,  the 
assumed  timber  harvesting  practice  also 
affected  evapotranspiration  in  two  respects. 
First,  during  the  snow  accumulation  and  melt 
seasons,  evaporation  from  the  snowpack  in 
the  small  openings  was  higher,  resulting  in 
greater  moisture  losses  than  from  snow  in 
uncut  forest.  Secondly,  evapotranspiration 
and  interception  losses  were  reduced  in 
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proportion  to  the  amount  of  forest  cover 
removed.  This  reduced  evapotranspiration 
resulted  in  lower  soil  water  deficits  on  the 
basin.  The  net  effect  was  an  overall  reduction 
in  evapotranspiration  and  resultant  in- 
creased water  yields.  Simulated  data  for  1947- 
71  water  years  are  shown  in  table  1.  (Note 
that  with  the  exception  of  snowpack  water 
equivalent,  all  hydrologic  components  are 
plotted  as  6-day  means  in  fig.  6.) 

The  simulated  average  runoff  increase 
for  the  1947-71  record  period  was  2.2  inches, 
which  resulted  from  a  2.2-inch  decrease  in 
evapotranspiration  losses,  with  no  change  in 
storage  during  the  average  water  year. 
Average  soil  water  requirements  on 
September  30  were  decreased  by  1.1  inches. 
As  discussed  above,  snowmelt  timing  and 
resultant  streamflow  were  also  changed. 
From  figure  6b,  it  is  seen  that  generated 
runoff  was  increased  during  April,  May,  and 
the  first  part  of  June  and  diminished 
somewhat  thereafter.  Because  the  generated 
flows  in  figure  6c  are  routed  through  natural 
storage  in  the  watershed  to  produce  the 


hydrograph,  it  is  reasonable  to  expect  that 
the  recession  limb  of  the  seasonal 
hydrograph  would  not  be  significantly 
changed  due  to  treatment.  However,  on  the 
rising  limb,  stream  discharges  would  be 
higher,  as  observed  from  watershed  studies  in 
the  area. 


Table  1 . --S imul ated  hydrologic  changes  resulting 
from  timber  harvesting  on  Deadhorse  Creek, 
Fraser  Experimental   Forest   (average  of  \3^7- 
71  water  years) 


Hydrologic 
component 


Water  balance 


Natural     Treated  Change 


Preci  p  i  tat  i  on 
Evapot ranspi  ration 
Soil  water  recharge 
requ i  rement 

beg  inning  (10/1) 

end  (9/30) 
Water  yield 


30.5 
16.8 


3.5 
3.5 
13.7 


Inahes  - 

30.5 
U.6 


0 

-2.2 


2.h 
15.9 


-1.1 
-1  .  1 
+2.2 


20 


o 
c 


-  15 


o 

3 


$  10 


O 

"o 
.0 


•  Before  treatment  ( 1943-55) 
A  After  treatment  (1956-72)  ^ 

A 

A 

A 

.  A 


April  -  Sept. 


Runoff 

Increase 

Since 

Harvest 

Year 

Inches 

1956 

3.3 

1957 

5.2 

1958 

3.7 

1959 

4.6 

1960 

5.4 

1961 

3.4 

1962 

3.9 

1963 

2.9 

1964 

3.5 

1965 

4.0 

1966 

2.6 

1967 

3.1 

1968 

1 .9 

1969 

2.5 

1970 

2.1 

1971 

3.6 

1972 

3.9 

Mean 

3.5  +0.8* 

*95  Percent  level  of 
confidence 


5  10  15  20 

E.St. Louis  Creek  runoff  (inches) 

Figure  5.  —  Pretreatment  and  posttreatment  correlations  of  seasonal  runoff  between  Fool  Creek 
and  adjacent  control  watershed,  Fraser  Experimental  Forest. 
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Long-Term  Hydrologic  Impacts  of 
Timber  Harvesting 

The  model  described  above  has  also  been 
used  to  simulate  the  long-term  effects  of  forest 
and  watershed  management  on  a  2,461-acre 
tributary  of  the  South  Tongue  River  in 
northcentral  Wyoming.  Pertinent  hydrologic 
characteristics  of  a  typical  forested  water- 
shed are  as  follows: 
Average  maximum  snowpack 

water  equivalent  15.5  inches 

Average  annual  precipitation      29.6  inches 


Average  annual  evapotran- 

spiration  15.8  inches 

Average  annual  runoff  13.8  inches 

Elevations  vary  from  8,000  feet  msl  to 
8,900  feet  msl.  Soils  are  derived  from  granitic 
rocks;  virtually  all  of  the  forest  cover  is 
lodgepole  pine.  To  illustrate  how  the  model 
was  used,  results  from  the  analysis  of  one 
planning  unit  will  be  summarized. 

In  addition  to  improving  water  yield,  the 
management  strategy  selected  for  this  exam- 
ple essentially  has  followed  recommen- 
dations published  by  Alexander  (1972),  which 
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Figure  6.  —  Simulated  average  water  balance  for  the  1947-71  water  years,  showing  changes 
resulting  from  patch-cutting  in  mature  subalpine  forest. 
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are  keyed  to  stand  descriptions,  insect  and 
disease  problems,  and  windfall  risk 
situations. 

Under  this  strategy,  all  of  the  old-growth 
timber  would  be  harvested  in  a  series  of 
patch-cuts  spread  over  a  planning  interval  of 
120  years.  At  intervals  of  30  years,  ap- 
proximately one-third  of  the  area  would  be 
harvested  in  small  openings  —  five  to  eight 
times  tree  height  —  distributed  over  the 
watershed  (table  2).  Forest  openings  would  be 
constructed  in  a  balanced  and  unified 
pattern  which  complements  the  natural 
landscape. 

Table  2 . --Watershed  management  strategy.  South 
Tongue  River  planning  unit,  Bighorn  National 
Forest,  Wyoming 


Management 
St  rategy 


Response  unit-' 


3 


Treatment 
Patch2 

Treatment 
Patch^ 
Thin  3 

Treatment 
Patch2 
Thin^ 

Treatment 
Thin  3 


1st  yr. 

31st  yr. 

61st  yr. 

91st  yr. 


^Each  unit  includes  1/3  of  the  total  area. 

^33  percent  of  area  occupied  by  openings  which 

are  5  to  8  tree-heights  in  diameter. 

^Thin  to  C.  ,, 
dmx/4 

As  each  one-third  of  the  old-growth  forest 
cover  on  the  planning  unit  is  patch-cut,  the 
forest  cover  density  on  the  previously  cutover 
areas  would  be  reduced  to  one-fourth  of  the 
natural  old-growth  forest  cover  density 
(Cdmx)-  At  the  end  of  the  planning  interval, 
all  of  the  openings  will  have  regenerated  and 
the  watershed  would  contain  groups  of  trees 
in  several  age  classes  from  reproduction  to 
those  ready  for  harvesting  on  the  originally 
cutover  areas.  The  management  strategy 
would  maintain  a  forest  cover  throughout  the 
planning  interval,  and  would  insure  suf- 
ficient seed  for  regeneration  from  trees  cut  on 
the  area,  or  standing  around  the  perimeters 
of  the  forest  openings. 

Projected  average  annual  water  yield 
increases  in  10-year  increments  under  this 
management  strategy  for  the  120-year  plan- 
ning interval  are  tabulated  in  table  3.  The 
increases  above  the  heavy  diagonal  line  in 
table  3  at  any  given  time  represent  the  overall 


response  resulting  from  preceding  manage- 
ment decisions.  The  data  below  the  line 
reflect  the  singular  effect  of  the  initial  patch- 
cut  on  one-third  of  the  planning  unit,  assum- 
ing that  it  were  the  final  decision  in  the 
strategy. 


Table  3. --Projected  changes  in  annual  water 
yield  resulting  from  timber  harvesting.  South 
Tongue  River  planning  unit,  Bighorn  National 
Forest,  Wyoming 


Interval 
(years) 


Water  yi/eld  increase,  by  treatment 


I  I  I 


IV 


0-10 

1 

33 

11 -20 

1 

59 

21-30 

0 

95 

31-40 

0 

74 

41-50 

0 

61 

51-60 

0 

51 

61-70 

0 

33 

71-80 

0 

08 

81-90 

0 

04 

91-100 

-0 

03 

101-1 10 

1 1 1-120 

Inches 


2.08 
2.07 
1  .72 


3.37 
2.86 
2.21 


2.90 
2.29 
1.76 


Water  yields  are  improved  throughout 
the  planning  interval,  with  the  highest  in- 
crease occurring  after  Treatment  III.  Pro- 
jected runoff  increases  in  relation  to  the 
pretreatment  base  period  during  each  treat- 
ment interval  are  as  follows: 


Treatment 


Runoff  Increase 
Percent 


I 

II 
III 

IV 


9.2 
14.2 
20.2 
16.8 


As  seen  from  table  3,  the  effect  of  the 
initial  patch-cut  (Treatment  I)  apparently 
persists  for  at  least  50,  and  perhaps  60  or 
more  years.  Thereafter,  the  effect  on  water 
yield  would  for  all  practical  purposes  be 
negligible.  (On  Fool  Creek,  in  central 
Colorado,  runoff  increases  have  not 
diminished  significantly  more  than  16  years 
after  strip  cutting.) 

The  projected  effects  of  timber 
harvesting  on  the  distribution  of  water 
available  for  streamflow  are  summarized  in 


11 


table  4.  These  values  represent  increments  of 
generated  runoff  and  not  routed  streamflow. 
Hence,  the  effects  of  watershed  storage  must 
be  considered  in  interpreting  the  data.  As 
seen  in  table  4,  inputs  from  snowmelt  are 
significantly  increased  during  April  and 
May,  and  diminished  in  June.  Minor  inputs 
to  streamflow  apparently  also  occur  in  July, 
while  none  occur  in  the  natural  state  due  to 
the  less  favorable  hydrologic  condition  of  the 
watershed. 


Table  4 . --Projected  changes  in  distribution  of 
water  available  for  streamflow,  South  Tongue 
River,  Bighorn  National  Forest,  Wyoming 


Average  change  in  run- 


III  IV 


Inches 


Apri 1  0.1  +0.9  +1.2  +1.2  +0.7 

May  7.5  +2.0  +2.2  +2.1  +1.3 

June  6.2  -2.8  -2.3  -2.7  -1.3 

July  0  +0.03  +0.05 


The  hydrologic  analysis  in  this  example 
indicated  that  the  magnitude  of  peak  flows 
would  be  changed  little  if  at  all  under  the 
proposed  management  strategy.  However, 
seasonal  peaks  would  occur  approximately 
one  week  earlier: 

Change  in 
peak  7-day 

generated         Change  in 

Treatment  runoff  timing  

 Inches  Days  


I  -0.3  -  9 

II  -0.5  -10 

III  -0.3  -  7 

IV  +0.3  -  5 


To  sum  up,  the  projected  overall 
hydrologic  impact  of  the  proposed  manage- 
ment strategy  would  be  to  increase 
streamflow  in  April  and  May  each  year 
throughout  the  120-year  planning  interval. 
This  accelerated  input  would  enlarge  early 
spring  flows  and  cause  the  hydrograph  to 
peak  approximately  one  week  earlier  than 
under  natural  conditions.  Hydrograph  peaks 
would  apparently  not  be  increased,  however, 
and  runoff  on  the  recession  side  of  the 
hydrograph  during  the  summer  months 
would  be  slightly  diminished. 


Discussion 

The  hydrologic  impacts  of  the  watershed 
management  practices  discussed  above  are 
but  two  examples  of  numerous  alternatives 
which  have  been  simulated  with  the  model 
described  in  this  paper.  The  model  has  been 
tested  and  calibrated  on  several  represen- 
tative drainage  basins  in  Colorado  and 
Wyoming.  The  areas  include: 

Wyoming:  South  Tongue  River,  Bighorn 
National  Forest;  East  Fork  of  the  Encamp- 
ment River,  Medicine  Bow  National  Forest. 

Colorado:  Eraser  River,  Arapahoe 
National  Forest;  Wolf  Creek,  San  Juan 
National  Forest. 


Timber  Models 

In  addition  to  improving  the  water  yields 
in  lodgepole  pine  and  spruce-fir  forests,  it 
should  be  emphasized  that  the  strategies 
selected  for  water  production  are  compatible 
with  the  conversion  of  old-growth  to  stands 
managed  from  the  regeneration  period  to 
final  harvest  for  timber  production.  Yield 
tables  that  report  probable  yields  of  wood 
that  result  from  specified  combinations  of 
site  quality,  frequency  and  intensity  of  thin- 
ning and  utilization  standards  provide  goals 
toward  which  conversion  can  be  directed. 
Procedures  for  deriving  yield  tables  for 
managed  stands  and  descriptions  of  the  main 
program  and  subroutines  have  been 
presented  for  lodgepole  pine  by  Myers  et  al. 
(1971)  (Program  LPMIST)  and  for  spruce-fir 
by  Alexander  et  al.  (1975)  (Program 
SPRYLD).  These  were  adapted  from  field 
and  computer  procedures  for  managed  stand 
yield  tables  originally  developed  by  Myers 
(1971). 

These  computer  programs  have  the 
capacity  of  producing  a  series  of  yield  tables 
which  show  how  projected  outcomes  will  vary 
in  response  to  changes  in  cultural  treatments 
and/or  variations  in  original  stand  and  site 
conditions.  Large  numbers  of  tables  each 
based  on  a  specific  set  of  alternatives  can  be 
computed  and  printed  at  the  cost  of  a  few 
cents  each.  This  provides  the  manager  with 
the  opportunity  to  examine  the  probable 
results  of  his  operations,  make  necessary 
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changes  in  management  goals,  and  study  the 
effect  of  these  changes  before  money  is  spent 
on  them  (Myers  1971). 


Linkage  Between  Hydrologic  and 
Timber  Models 

The  hydrologic  model  and  the  timber 
models  (LPMIST  and  SPRYLD)  are  linked 
by  means  of  the  forest  cover  density  variable 
(C^)  as  defined  previously.  Forest  cover  den- 
sity is  assumed  to  vary  as  a  function  of 
time  according  to  the  expression  (eq.  [8]): 

where 

=  intermediate  forest  cover  den- 
sity after  cutting  is  sufficient  to 
reestablish  the  stand.  When 
t  <  t^^,       =  0., 

<I>  =  the  time  in  years  from  tc2  at 
which  maximum  forest  cover 
density  is  reached.  This  param- 
eter will  vary  according  to  tree 
species,  local  environment,  and 
stand  condition,  and, 

^dmx  ~  maximum  (natural  old  growth) 
forest  cover  density  expressed  as 
a  decimal. 

In  the  hydrologic  model,  logging  cor- 
responds to  a  reduction  of  the  forest  cover 
density  (C{j).  Thus,  the  degree  that  C(j  is 
reduced  depends  on  the  relative  changes  in 
basal  area.  As  stated  previously,  in  old- 
growth  stands,  if  Cd  is  reduced  by  50  percent 
or  less  from  Cdmx>  it  is  assumed  that  forest 
cover  density  does  not  increase  subsequent  to 
cutting.  However,  if  is  reduced  more  than 
50  percent  from  Cdmx,  but  not  clearcut, 
equation  [8]  is  used  to  simulate  redevelop- 
ment of  the  canopy  with  time.  In  the  event 
that  Cd  is  reduced  to  zero  (clearcut),  is 
replaced  by  a  new  parameter,  in  equation 
[8],  which  then  computes  redevelopment  of 
the  canopy  with  time  under  "managed 
stand"  conditions. 

No  relationship  has  been  established 
between  Cd  and  basal  area,  average 
diameter,  and  site  index.  However,  calibra- 


tion studies  indicate  it  is  reasonable  to 
assume  that  a  given  percentage  reduction  in 
cover  density  corresponds  reasonably  well  to 
a  similar  reduction  in  basal  area.  Com- 
parisons between  basal  area  levels  after 
patch-cutting  and  forest  cover  density  are 
summarized  for  a  specified  set  of  stand  con- 
ditions in  table  5. 

Table  5 . ""Compa r i sons  of  basal  area^  after 

initial  clearcutting  with  forest  cover  density 


Basal  area  Forest  cover 

Years  since  (ft^/ac)  density  (C^) 

initial  cut     

Before    After  Before  After 


Lodge pole 

p  i  ne  ( 

(})'  =  30  years) 

0 

130 

0 

0.30  0 

30 

no 

61 

.30  0.16 

60 

129 

92 

.30  .21 

90 

100 

.30  .21 

120 

]kG 

0 

.30  0 

Spruce- 

fir 

=  60  years) 

0 

325 

0 

.55  0 

50 

9't 

63 

.38  .2k 

80 

153 

99 

.55  .35 

110 

176 

100 

.55  .31 

140 

160 

0 

.55  0 

Basal  areas  for  lodgepole  pine  computed  from 
Program  LPMIST  for  site  index  60,  and  initial 
and  subsequent  growing  stock  levels  of  100. 


Basal  areas  for  spruce-fir  computed  from 
Program  SPRYLD  for  site  index  80,  and  initial 
and  subsequent  growing  stock  levels  of  100. 


Applications 

Lodgepole  Pine.  —  The  management 
strategy  (table  6)  is  similar  to  that  in  table  2 
for  the  Bighorn  National  Forest,  but  for 
another  area  in  central  Colorado. 

Similar  to  the  watershed  management 
strategy  in  the  previous  example,  water 
available  for  streamflow  on  each  response 
unit  is  substantially  increased  by  the  patch- 
cutting  of  lodgepole  pine  (table  7,  fig.  7).  In  the 
cleared  areas,  the  increase  in  water  available 
for  streamflow  is  maintained  at  a  higher  level 
by  frequent  intermediate  thinnings 
throughout  the  rotation  than  if  the  forest  was 
allowed  to  return  to  preharvest  conditions 
naturally.  Figure  8  shows  projected  water 
yield  changes  from  initial  patch-cutting  and 
subsequent  thinning  on  one  response  unit.  In 
the  example  used  here,  the  forest  manager 
also  wishes  to  determine  the  growing  stock 
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levels  that  will  maximize  volume  production 
in  board  feet  within  the  limits  imposed  by  the 
cutting  strategy  for  water  production.  Since 
water  yield  is  unaffected  by  site  quality,  an 
average  site  index  (60)  has  been  chosen. 
Furthermore,  simulation  analyses  indicate 
that  water  yields  are  little  affected  by  any 
combination  of  initial  and  subsequent  grow- 
ing stock  levels  in  managed  stands  that 
range  from  <  80  to  >  120  ftVac.  Length  of 
rotation  is  120  years  with  a  30-year  cutting 
cycle.  Alternatives  that  call  for  more  than  one 
precommercial  thinning  are  unacceptable. 
Minimum  commercial  volumes  per  acre  are 
400  cubic  feet  and  1500  board  feet.  The 
manager  expects  that  his  procedures  for 
regenerating  each  patch-cut  will  result  in  a 
new  stand  that  contains  1,000  trees  per  acre 
by  age  30,  with  an  average  stand  diameter  of 
4.5  inches.  Furthermore,  dwarf  mistletoe  in- 
fection will  not  occur  during  the  life  of  the 
stand. 

Table  6 . --Wate rshed  management  strategy  for 
lodgepole  pine,  Deadhorse  Creek  planning 
unit,  Fraser  Experimental  Forest,  Colorado^ 


Lodgepole  pine 


Management 
strategy 


Response  unit 


Treatment 

Patch2 
Treatment 

Patch2 

Thin  3 
Treatment 

Patch2 

Thin3 
Treatment 

Thin3 
Treatment 

Harvest^ 

Thin^ 


I  I  I 

IV 

V 


1st  y r. 

31st  yr. 

61st  yr. 

91st  yr. 

121st  yr. 


^Aspect:  SSE 

Elevation:     10,500  ft.  msl 

Slope:     30  percent 
^33  percent  of  the  area  cut  in  openings  5  to  8 

times  tree  height. 

^Thin  to  growing  stock  levels  specified  under 
section  on  Timber  Management  Alternatives. 

A  few  of  the  yield  tables  produced  by 
LPMIST  are  reproduced  in  Appendix  A.  For 
the  situation  described  above,  yields  and 
number  of  precommercial  thinnings  are  of 
greatest  interest.  These  items  are  sum- 
marized in  tables  8  and  9  for  the  9  yield  tables 
produced.  Only  the  combination  of  low  initial 
and  low  subsequent  growing  stock  levels 
meets  the  requirement  of  only  one 
precommercial  thinning. 


0  20         40         60         80         100        120  140 

Time  since  initial  cutting  (years) 

Figure  7.  —  Projected  water  yield  changes 
from  management  strategies  outlined  in 
tables  7  and  11. 


Table  7. ""Projected  changes  in  annual  water 
yield  resulting  from  tiber  harvesting  in 
lodgepole  pine,  Deadhorse  Creek  planning 
unit,  Fraser  Experimental  Forest,  Colorado 


Interval 
(years) 


Water  yield  increase,  by  treatment 


III  IV 


0-10 
1  1"20 
21"30 

41-50 
51"60 

61"70 
71"80 
8l"90 

91"100 
101-110 
1 11-120 


Inches 


1.72 
1.73 
1 .29 


l.kk 
2.30 
\.5h 


3.12 
2.86 
2.11 


2.8it 
1 .98 
1 .38 


Additional  comparisons  can  be  made  to 
include  such  factors  as  probable  thinning 
costs,  cubic  yields  from  thinnings  not  com- 
mercial for  board  feet,  and  the  average  size  of 
tree  produced.  As  expected,  the  current  crop 
produces  more  board  feet  in  120  years  at  high 
subsequent  levels  of  growing  stock,  but  two 
precommercial  thinnings  would  be  required. 

Spruce-fir.  —  Water  available  for 
streamflow  on  each  response  unit  is  also 
substantially  increased  by  patch-cutting 
spruce-fir  in  small  openings  5  to  8  times  tree 
height  according  to  the  management 
strategy  outlined  in  table  10  (table  11,  fig.  7). 
However,    on    the    cleared   areas,  water 
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available  for  streamflow  decreases  after  the 
initial  cutting  at  about  the  same  rate  whether 
or  not  intermediate  thinnings  are  made 
throughout  the  rotation.  This  is  illustrated  by 
figure  9,  which  shows  projected  water  yield 
changes  from  initial  patch-cutting  and  subse- 
quent thinning  on  one  response  unit  com- 
pared with  no  thinning  after  initial  harvest. 
There  are,  however,  other  advantages  to 
thinning  spruce-fir.  Growth  is  concentrated 
on  fewer  stems,  and  total  yields  of  usable 
products  are  increased.  In  the  example  here, 
the  forest  manager  also  wishes  to  determine 
the  growing  stock  levels  that  will  maximize 
volume  production  in  board  feet  within  the 
limits  imposed  by  the  watershed  manage- 
ment strategy.  Water  yields  are  little  affected 
in  spruce-fir  forests  by  either  site  quality  or 
the  growing  stock  levels  that  are  likely  to  be 
timber  management  goals.  An  average  site 
index  of  80  was  chosen.  Length  of  rotation  is 
120  years  (breast  height  age)  with  a  30-year 
cutting  cycle. 2  Alternatives  that  call  for  more 
than  one  precommercial  thinning  are  un- 
acceptable. Minimum  commercial  volumes 
per  acre  are  400  cubic  feet  and  2,000  board 
feet.  The  manager  expects  that  his  procedure 
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~  Lodgepole  pine 
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 Unmanaged 
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20         40         60         80  100 
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Figure  8.  —  Projected  water  yield  changes  for 
managed  (all  options)  and  unmanaged 
lodgepole  pine  following  initial  patch- 
cutting  on  one  response  unit.  Cutting  cycle 
30  years;  site  index  60  feet;  all  initial  and 
subsequent  stocking  levels  likely  to  be 
timber  management  goals. 


2  Age  in  SPRYLD  is  age  at  breast  height.  The  50-year 
interval  between  patch-cutting  and  the  first  thinning  in 
table  11  is  to  allow  a  minimum  of  20  years  for  spruce  and 
fir  trees  to  regenerate  and  grow  to  4.5  feet  in  height.  A 
120-year  rotation  is  therefore  at  least  140  years  in  the 
total  age  of  the  stand. 


for  regenerating  each  patch-cut  will  result  in 
a  new  stand  that  contains  850  trees  per  acre 
with  an  average  stand  diameter  of  4.5  inches 
by  b.h.  age  30  years. 


Table  8. --Number  of  precommercial  tliinnings 
based  on  minimum  board  feet  volumes,   if  each 
of  ttie  9  combinations  of  initial  and  subse- 
quent growing  stocl<  levels  is  executed  as 
specified  by  the  data  decks  for  spruce-fir 
(SPRYLD)  and  lodgepole  pine  (LPMIST) 


Initial  thinning 
basal  area  level 
(ftVac) 


Subsequent  basal  area 
level  (ft2/ac) 


80 


100 


120 


Spruce-fir  (Program  SPRYLD) 


80 
100 
120 


Lodgepole  pine  (Program  LPMIST) 


80 
100 
120 


Table  9. — Yields  in  thousand  board  feet,  in- 
cluding commercial  thinning  of  9  combinations 
of  initial  and  subsequent  growing  stock  levels, 
spruce-fir  and  lodgepole  pine 


Initial  th  i  nn  i  ng 
basal  area  level 

Subsequent  basal 
level  (ft^/ac) 

a  rea 

(ft  Vac) 

80  100 

120 

Spruce-f i 

r  (Program  SPRYLD) 

80 

32.7  36.2 

'»1.3 

100 

33.3  37.'* 

/to. 6 

120 

33.5  37.6 

hO.k 

Lodgepole  p 

ine  (Program  LPMIST) 

80 

2't.2  25.0 

27.5 

100 

21.0       2'). 8 

27.3 

120 

20.8  2'4.0 

27.0 

A  few  of  the  yield  tables  produced  by 
SPRYLD  are  reproduced  in  Appendix  B.  For 
the  situation  described  above,  yields  and 
number  of  precommercial  thinnings  are  sum- 
marized in  tables  8  and  9  for  the  9  tables 
produced.  The  combination  of  high  initial 
and  intermediate  subsequent  growing  stock 
levels  produces  the  greatest  volume  with  one 
precommercial  thinning. 

Comparisons  can  also  be  made  to  include 
thinning  costs,  cubic  yields  from  thinnings 
not  commercial  for  board  feet  and  the  average 
size  of  tree  produced.  The  current  crop 
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produces  more  board  feet  with  the  combina- 
tion of  low  initial  and  high  subsequent  grow- 
ing stock  levels,  but  two  precommercial  thin- 
nings will  be  required. 


■§2 


-  0 


Spruce-fir 


Managed 
Unmanaged 


20 


40         60         80  100 
Time  since  initial  cutting  (years) 


120 


140 


Figure  9.  Projected  water  yield  changes  for 
managed  (all  options)  and  unmanaged 
spruce- fir  following  initial  patch-cutting  on 
one  response  unit.  Cutting  cycle  30  years 
(50  years  initial  cut);  site  index  80  feet;  all 
initial  and  subsequent  stocking  levels  like- 
ly to  be  timber  management  goals. 


Table  1 1 . --Projected  changes  in  water  yield 
resulting  from  timber  harvesting  in  spruce- 
fir,  Deadhorse  Creel<.  planning  unit,  Fraser 
Experimental  Forest,  Colorado 


Interval 
(years) 


Water  yield  increase,  by  treatment 


Inches 


0-10 

0.91 

1  I  -20 

1  . 20 

21  -30 

31-40 

.10 

41-50 

.10 

51-60 

1.26 

61-70 

1 .49 

71-80 

1.21 

81-90 

0.71 

91-100 

.42 

101-110 

1  .22 

1 1 1-120 

2 

121-130 

1 

131-140 

I 

.54 
.09 


Table  1 0 . --Watershed  management  strategy  for 
spruce-fir,  Deadhorse  Creek  planning  unit, 
Fraser  Experimental  Forest,  Colorado^ 


strategy 


Treatment 

1 

Patch2 

1st 

yr. 

X 

Treatment 

1  1 

Patch^ 

51st 

yr. 

X 

Thin  3 

X 

Treatment 

1  1  1 

Thin  3 

81st 

yr. 

X 

Treatment 

IV 

Patch^ 

101st 

yr. 

Thin  3 

X 

Treatment 

V 

Thin  3 

1 1 1th 

yr. 

X 

Treatment 

VI 

Thin  3 

131st 

yr. 

X 

Treatment 

V  1  1 

Harvest^ 

141st 

yr. 

X 

^Aspect : 
E 1 evat  i  on 

NE 
: 

0,200 

ft.  msl 

Slope:     35  percent 

■33  percent  of  the  area  cut  in  openings  5  to  8 
times  tree  height. 

'Thin  to  growing  stock  levels  specified  under 
section  on  Timber  Management  Alternatives. 
Yield  tables  for  spruce-fir  are  based  on  breast 
height  age.     On  a  30-year  cutting  cycle  a 
minimum  of  20  additional  years  will  be  required 
for  trees  to  reach  breast  height. 


Summary 

The  following  paragraphs  highlight  the 
practical  aspects  of  this  work,  and  summarize 
important  principles  which  should  be  con- 
sidered in  land  use  planning. 

•  Highest  water  yields  result  when  old- 
growth  subalpine  forests  are  harvested  in 
small  patches.  When  forest  openings  are:  (1) 
less  than  8  tree-heights  in  diameter;  (2) 
protected  from  wind,  and  (3)  interspersed  so 
that  they  are  5  to  8  tree-heights  apart,  an 
optimum  pattern  of  snow  accumulation 
results.  More  snow  is  deposited  in  the 
openings,  and  less  snow  accumulates  in  the 
uncut  forest  so  that  total  snow  on  headwater 
basins  is  not  significantly  increased. 

•  Snowmelt  in  the  small  openings  on  all 
aspects  is  more  rapid  than  in  the  uncut  forest. 
This  accelerated  melt  causes  streamflow  to  be 
higher  on  the  rising  limb  of  the  hydrograph 
than  before  harvest  cutting.  When  there  is 
considerable  natural  regulation  in  the  form  of 
deep  porous  soils,  recession  flows  should  not 
be  changed  appreciably  and  annual  flood 
peaks  are  not  significantly  increased  provid- 
ed that  the  forest  cover  on  no  more  than  50 
percent  of  the  watershed  is  removed  in  a 
system  of  small  openings. 
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•  Simulation  analyses  indicate  that  un- 
der a  patch-cut  alternative,  water  yield  in- 
creases on  south  slopes  are  at  least  as  large  as 
corresponding  increases  from  north  aspects. 
Hence,  there  is  no  reason  to  favor  areas  with 
the  highest  natural  water  yield  if  the  objec- 
tive is  to  maximize  water  yield  from  old- 
growth  subalpine  forests. 

•  Due  to  the  considerable  length  of  time 
that  it  takes  for  subalpine  coniferous  forests 
to  regenerate,  increased  water  yields  from 
patch-cutting  can  go  undiminished  for  30 
years  and  longer.  Even  after  this  period  of 
time,  it  is  conceivable  that  30  additional  years 
will  be  required  before  runoff  increases  from 
the  initial  timber  harvest  are  completely 
erased. 

•  It  should  be  emphasized  that  the 
pattern  in  which  trees  are  harvested  deter- 
mines whether  or  not  runoff  will  be  increased. 
For  example,  when  the  forest  cover  is  remov- 
ed in  large  clearcut  blocks  or  by  selective 
cutting  of  individual  trees,  increased  water 
yields  will  be  far  less  than  that  attained  if  the 
same  volume  of  timber  is  harvested  in  a 
system  of  small  dispersed  forest  openings. 
Under  some  conditions,  streamflow  may  ac- 
tually be  decreased  when  timber  is  selectively 
harvested  or  clearcut  in  large  blocks. 

•  In  much  of  the  Rocky  Mountain  Region, 
timber  harvesting  which  produces  the  most 
additional  water  is  ecologically  sound.  If 
done  properly,  it  does  not  reduce  water  quali- 
ty; it  is  a  silviculturally  acceptable  procedure 
and  compatible  with  the  guidelines  recently 
developed  from  research  in  old-growth  sub- 
alpine forests  (Alexander  1972,  1973,  1974). 
The  strategies  selected  for  optimum  water 
production  are  compatible  with  the  conver- 
sion of  old-growth  to  stands  managed  from 
the  regeneration  period  to  final  harvest  for 
timber  production. 

•  Procedures  are  available  for  projecting 
long-term  yields  of  both  wood  and  water 
resulting  from  a  broad  array  of  management 
alternatives.  This  capability  provides  the 
manager  with  the  opportunity  to  examine  the 
probable  results  of  his  operations  within  the 
context  of  multi-resource  management. 


Conclusion 

One  of  the  major  shortcomings  of  the 
models  described  in  this  paper  is  the  lack  of 


sufficient  data  for  validation,  particularly 
with  respect  to  man's  long-term  impacts  on 
the  timber  and  water  resources.  Nevertheless, 
the  models  produce  expected  results  based  on 
experience  and  the  state-of-the-art.  It  is 
believed  that  the  output  from  the  examples 
above  contain  the  type  of  information  which 
hydrologists,  silviculturalists,  and  land  use 
planners  need  to  know  in  order  to  make 
difficult  management  decisions.  The  ability 
of  the  models  described  in  this  paper  and 
other  similar  models  to  integrate  complex 
forest  and  water  systems  make  them  unique 
and  powerful  tools  for  evaluating  a  broad 
array  of  land  management  alternatives. 
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Appendix  A 
Typical  Yield  Tables  Produced  by  LPMIST 

YIELDS  PER  ACRE  OF  EVEN-AGED  STANDS  OF  LODGEPOLE  PINE 

SITE  INDEX  60 
THINNING  INTENSITY-  INITIAL-       80.     SUBSEQUENT-  80. 

ENTIRE  STAND  BEFORE   AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 

(YEARS ) 

TREES 
NO  . 

BASAL 
AREA 
SO .FT . 

AVERAGE 
D  .B.H. 

AVERAGE 
HEIGHT 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
C  U  •  F  T  • 

SAMTIMBER 
VOLUME 
BD.FT. 

TREES 
NO  . 

BASAL 
AREA 
SQ  .F  T  . 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAHTIMBER 
VOLUME 
BD.FT. 

30 
30 

10  0  0 
325 

110 

50 

^  .  5 
5.3 

2't 

26 

127  0 
630 

170 
170 

g 

0 

675 

60 

61.0 

0 

0 

«0 

322 

72 

6.1* 

32 

116  0 

720 

0 

50 

321 

93 

7.3 

36 

1670 

13i|0 

0 

60 
60 

321 
105 

115 

8.1 

8  •  7 

l<2 

U  3 

2'*20 
165  0 

2120 
li»90 

8700 
6100 

136 

39 

770 

630 

2600 

70 

182 

91 

9.6 

1)8 

2220 

2050 

8600 

80 

179 

106 

10. <i 

52 

2790 

2620 

11100 

90 
90 

179 
10<i 

122 
80 

11.2 
11.9 

56 
57 

3ii80 
2280 

3280 
2160 

11.200 
9600 

75 

1(2 

1200 

1120 

1.600 

100 

10<i 

93 

12.8 

60 

271.0 

2600 

11800 

IID 

10  . 

106 

13.7 

63 

32i»0 

3090 

1<«300 

120 

lO'- 

119 

111. 5 

65 

3730 

3560 

16800 

DMARF 

MISTLETOE 

INFECTION 

DID  NOT 

OCCUR  DURING  THE  ROTATION  OF 

120.  YEARS. 

CERCH. 

CU.  FT.  ■ 

-  TREES  6. 

0  INCHES 

O.B.H.  ANO 

LARGER  TO 

I.-INCH  TOP. 

BO.   FT.  -  TREES 

6.5  IKCHES  O.B.H. 

AND  LARGER 

TO  6-INCH 

TOP. 

YIELDS  PER  ACRE  OF  EVEN-AGED  STANDS  OF  LODGEPOLE  PINE 

SITE  INDEX  60 
THINNING  INTENSITY-  INITIAL-     120.     SUBSEQUENT-  120. 

ENTIRE  STAND  BEFORE  AND  AFTER  THINNING  PERIODIC   INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.e.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAHTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAHTIMBER 
VOLUME 

eo.FT. 

30 
30 

1000 
505 

110 

72 

ii.S 
5.1 

21. 

25 

1270 
880 

180 
ISO 

0 
0 

1.95 

38 

390 

0 

0 

1.0 

502 

99 

6.0 

32 

151.0 

760 

0 

50 

500 

122 

6.7 

35 

211.0 

151.0 

0 

60 
60 

1.98 
322 

11.5 
107 

7.3 
7.8 

1.1 

1.2 

2990 
2250 

21.10 
1920 

0 
0 

176 

38 

71.0 

^90 

0 

70 

322 

127 

8.5 

1.7 

3010 

2700 

11100 

80 

;22 

11.5 

9.1 

52 

3770 

31.70 

11.300 

90 
90 

^21 
'07 

165 
120 

9.7 
10.3 

55 
56 

1.580 
3370 

1.21.0 
3170 

17900 
13300 

111. 

1.5 

1210 

1070 

1.600 

100 

'07 

137 

11.0 

59 

1.070 

381.0 

16500 

110 

'07 

152 

11. E 

62 

1.670 

1.1.10 

191.00 

120 

'07 

168 

12.2 

61. 

5300 

5020 

22400 

DWARF 

MISTLETOE 

INFECTION 

DID  NOT 

OCCUR  DURING  THE  ROTATION  OF 

120.  YEARS. 

MERCH. 

CU.  FT.  - 

TREES  6. 

0  INCHES 

D.B.H.  AND 

LARGER  TO 

1.-INCH 

rop. 

BD.  FT 

.  -  TREES 

6.5  INCHES  D.B.H. 

ANO  LARGER 

TO  6-INCH 

TOP. 
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Appendix  B 
Typical  Yield  Tables  Produced  by  SPRYLD 


YIELDS  PER   «CRE  OF  HftNAGEOi   EVEN-AGED  STANDS  OF  EN6ELMANN  SPRUCE  AND  SUBALPINE  FIR 
SITE  INDEX     80,     30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS=  INITIAL  -       80.,   SUBSEQUENT  -  80. 
ENTIRE  STAND   EEFORE  AND  AFTER  THINNING  PERIODIC  INTERMEDIATE  CUTS 


STAND 
AGE 
(YEARS) 

TREES 
NO. 

BASAL 
ARE  A 
SQ.FT. 

AVERAGE 
0  .  B  •  H  . 
IN. 

AVERAGE 
HE IGHT 
FT. 

TOTAL 

V  U  LU  nL 

CU.FT. 

MERCHANT- 

AQI  C     11  m  IlUC 

AdLC  vULUnt 

CU.FT. 

SAHTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

30 
30 

850 
31<i 

52 

<t.5 
5.5 

28 
29 

1010 
680 

3<tO 
3i«0 

0 
0 

536 

1,2 

330 

0 

itO 

311 

78 

6.8 

38 

1260 

900 

0 

50 

305 

106 

8.0 

2030 

1650 

>i200 

60 
60 

295 
l<t9 

133 
80 

9.1 
9.9 

53 
53 

2920 
1310 

2530 
1610 

8300 
5  7  0  0 

5  3 

1110 

92  0 

70 

11(5 

99 

11.2 

59 

2510 

2310 

9200 

80 

l<t5 

122 

IZ.k 

65 

3370 

3150 

13600 

90 
90 

1'.5 

55 

1<|6 
80 

13.6 
16.3 

70 
71 

2<t<t  0 

<tl20 
231,0 

18700 
11100 

90 

66 

1910 

1780 

100 

55 

95 

17.8 

75 

300  0 

2910 

liiSOO 

110 

55 

112 

19.3 

79 

3610 

3520 

18300 

120 

55 

130 

20.8 

82 

'•290 

<i2  0  0 

TOTAL 

22500 
YIELDS 

76I|0 

6900 

MINIMUH 

CUTS 

FOR  INCLUSION  IK 

TOTAL   YIELDS — 

<i00.   CUBIC  FEET 

AND  2000. 

BOARD 

FEET 

HERCH. 

CU.   FT.  -  TREES 

5.0   INCHES  O.e.H. 

AND  LARGER  TO  i»-INCH 

TOP. 

SAHTIKEER 
VOLUME 
ED. FT. 


2600 


7600 


32700 


BD.   FT.  -  TREES  8.0   INCHES  D.B.H.   AND  LARGER  TO  6-INCH  TOP. 


YIELDS  PER  ACRE  OF   MANAGED,   EVEN-AGED  STANDS  OF  ENGELMANN  SPRUCE  AND  SUBALPINE  FIR 
SITE   INDEX     80,     30-YEAR  CUTTING  CYCLE 
THINNING  LEVELS^   INITIAL  -     120.,  SUBSEQUENT  -  120. 
ENTIRE  STAND  EEFORE  AND  AFTER  THINNING  PERIODIC  INTERMEDIATE  CUTS 


STAND 
AGE 
•YEARS) 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

AVERAGE 
D.B.H. 
IN. 

AVERAGE 
HEIGHT 
FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAWTIMBER 
VOLUME 
BD.FT. 

TREES 
NO. 

BASAL 
AREA 
SQ.FT. 

TOTAL 
VOLUME 
CU.FT. 

MERCHANT- 
ABLE VOLUME 
CU.FT. 

SAHTIMBER 
VOLUME 
BD.FT. 

30 
30 

850 
505 

72 

't.5 
5.1 

28 
29 

1010 

870 

31)0 
3<<0 

0 
0 

31.5 

22 

11.0 

0 

0 

dO 

'»98 

108 

6.3 

37 

1620 

1050 

0 

50 

If  85 

I'd 

7.3 

kS 

25',0 

191(0 

0 

60 
60 

<»66 
259 

175 
117 

8.3 
9.1 

52 
52 

3660 
2570 

30i<0 
2230 

8700 
7300 

207 

58 

1090 

810 

11.00 

70 

2it9 

l<tl 

10.2 

59 

3^80 

3130 

11800 

80 

2<»9 

173 

11.3 

i>680 

Ii300 

17800 

90 
90 

2*t9 
103 

205 
120 

12.3 
l<t.6 

69 
70 

5960 
3600 

5570 
31,30 

2'i'tOO 
15900 

11.6 

85 

2360 

211.0 

8500 

100 

103 

l^ 

15.9 

75 

1.360 

20900 

110 

103 

166 

17.2 

78 

SklO 

5220 

26300 

120 

103 

190 

18. <i 

82 

6310 

6120 

TOTAL 

31900 
YIELDS 

9900 

9070 

1.01.00 

MINIMUM 

CUTS 

FOR  INCLUSION  IN 

TOTAL  YIELDS—     <.00.   CUBIC  FEET 

AND  2000. 

BOARD 

FEET 

KERCH. 

CU.    FT.  -  TREES  5.0   INCHES  D.E. 

H.  AND  LARGER  TO  <»-INCH 

TOP. 

BD.  FT. 

-  TREES  8.0 

INCHES  D.B.H.  AND 

LARGER  TO 

6-INCH  TOP. 
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